Device-independent demonstration of genuine multipartite entanglement 
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Entanglement in a quantum system can be demonstrated experimentally by performing the mea- 
surements prescribed by an appropriate entanglement witness [T]. However, the unavoidable mis- 
match between the implementation of measurements in practical devices and their precise theoretical 
modelling generally results in the undesired possibility of false-positive entanglement detection [2]. 
Such scenarios can be avoided by using the recently developed device-independent entanglement 
witnesses (DIEWs) for genuine multipartite entanglement [3]. Similarly to Bell inequalities, DIEWs 
only assume that consistent measurements are performed locally on each subsystem. No precise 
description of the measurement devices is required. Here we report an experimental test of DIEWs 
on up to six entangled 40 Ca + ions. We also demonstrate genuine multipartite quantum nonlocality 
between up to six parties with the detection loophole closed. I 



Entanglement enables many quantum tasks, includ- 
ing scalable quantum communication [4j , secure quantum 
key distribution [5J [5] and quantum computing [7] . En- 
tanglement thus represents a key resource for quantum 
information processing. As the experimental systems de- 
veloped to host and manipulate entangled states diver- 
sify [81-11 1 j . clear tests are desired that are able to certify 
properly the entanglement of such states. In particular, 
tests of genuine multipartite entanglement are needed to 
show that the entanglement truly involves all of their 
constituents PQ. 

Tools to detect genuine multipartite entanglement have 
been developed pQ and used to demonstrate, for instance, 
entanglement between up to 8-10 photonic qubits [T2IIT5] 
or 14 ions [14] . Namely, it was shown that every multipar- 
tite entangled state can be detected by a hermitian oper- 
ator known as an entanglement witness, whose expecta- 
tion value is negative on the considered state but positive 
on all biseparable states. Any measurement of a witness 
yielding a negative value thus certifies that the measured 
state is multipartite entangled. Experimentally speak- 
ing, entanglement witnesses are convenient because they 
can be expressed in terms of simple measurable observ- 
ables specific to the physical system, in contrast to joint 
measurements over multiple copies. 

Although entanglement witnesses are very well under- 
stood from a theoretical point of view, their applica- 
tion to practical systems strongly relies on the way that 
concrete systems are modeled. For example, entangle- 
ment witnesses generally take different forms depending 
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on the dimension of the Hilbert space under considera- 
tion. An artificial truncation of this space can thus lead 
to wrong entanglement detection (e.g., in photon num- 
ber |15j). Similarly, any mismatch between the measure- 
ments prescribed and those performed, such as system- 
atic errors, can also lead to flawed entanglement detec- 
tion [2]. Entanglement witnesses can also be evaluated 
on a state reconstructed by tomography [TBI [17] , but the 
reconstruction also relies on a detailed model of the mea- 
surement devices and Hilbert space dimensionality (see 
e.g. Ref. [H] for possible consequences of using a wrong 
model). In particular, quantum state tomography is also 
sensitive to systematic measurement errors [2]. 

Systematic errors and untenable assumptions may lead 
one to conclude that any proper experimental entan- 
glement certification requires a particularly good prior 
knowledge of the measurement devices. But this is not 
the case: a rudimentary measurement model, whose va- 
lidity can in principle be checked in practice, is sufficient 
to demonstrate entanglement properties of quantum sys- 
tems. Namely, entanglement can be certified under the 
sole hypothesis that each subsystem can be addressed 
individually with independent measurement devices. In- 
deed, this condition allows one to ensure that, upon mea- 
surement of a separable quantum state, the raw statis- 
tics satisfy all Bell inequalities pQ. Violation of a Bell 
inequality under these conditions thus certifies that the 
measured state was entangled. In contrast with conven- 
tional entanglement witnesses, this conclusion holds in- 
dependently of the internal mechanisms of the measure- 
ment devices used in the experiment [TH] (device indepen- 
dently as in recent tests of classical and quantum dimen- 
sions (2HH23)- Because Bell inequalities, in general, do 
not detect genuine multipartite entanglement, dedicated 
device- independent witnesses were recently developed in 
the form of DIEWs [1 [23 EH] • 
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Several architectures have realized multi-particle en- 
tanglement [5HTT] but trapped ions have demonstrated 
particularly high fidelity in the preparation, manipula- 
tion and detection of quantum states. Similar to other 
systems, however, such ion trap architectures are affected 
by crosstalk. This imperfection occurs on systems in 
which measurements or operations intended on individ- 
ual qubits affect neighboring qubits (close physically or 
in frequency, for example). Because crosstalks were not 
considered in the original derivation of DIEWs [3J, we 
adapted the bound of the tested witnesses to take this 
imperfection into account. We then demonstrate gen- 
uine multipartite entanglement between three, four and 
six trapped ions with the aid of these DIEWs. 

As described in detail in Ref . [3J , let us consider an en- 
tangled state consisting of n parties or subsystems and a 
device-independent scenario where the measurement set- 
tings and outcomes are only refered to by indices and 
without alluding to the specifics of how these measure- 
ments are performed. We consider that a party j of this 
state is measured in one of m possible ways. We thus 
index this measurement setting by Sj = 0, . . . , m — 1. 
We also consider that the local measurement of each 
party yields one of two outcomes that we index by 
Tj = 0,1. We can thus identify a measurement on the 
entire n-partite system by the vector s — (s\, . . . , s n ) £ 
{0, . . . ,171— 1}". Similarly, we describe each possible out- 
come of a measurement on the entire system by a vector 
r = (ri, . . . , r n ) £ {0, 1}™. After measuring the state, the 
conditional probability, or raw statistics, of observing the 
outcome r given that the choice of measurement settings 
s is applied to each of the n qubits is denoted by P(r|s). 
For an n-partite state with the described measurements 
and outcomes, we define the DIEW parameter as 

s=0 (mod m) 

£ (-l)W m Eg (1) 

s=l (mod m) 

where Eg = l) r P(r\s) is the n-partite correlator; 

s = J2j s j and r = £V r j ( see R- ei - Hi)- 



P(r\s) = Tr 



As shown in Ref. 
correlations satify 



because all biseparable quantum 
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this inequality can be used as a DIEW to detect gen- 
uine n-partite entanglement. In contrast with usual en- 
tanglement witnesses, I nm provides no explicit descrip- 
tion of the observables to be measured on each sub- 
system. Yet, every practical evaluation of Eq. ([!]) is 
necessarily performed with some measurement operators 
Ml i so that the observed correlations can be written as 



Here M 3 , 



is the measure- 



ment operator acting locally on subsystem j correspond- 
ing to setting Sj and outcome rj, and p is the measured 
quantum state. The DIEW parameter can thus be un- 
derstood as the expectation value of a quantum operator 
T nm whose biseparable bound is given by B nm (see Sup- 
plementary Information for a definition of T nm ) . In par- 
ticular, the quantity W nm = B nm l — I nm is a standard 
entanglement witness for genuine multipartite entangle- 
ment pQ . Therefore, any practical evaluation of a DIEW 
amounts to test a standard witness that is adapted to the 
measurement operators available during the experiment. 
In addition, DIEWs need less measurements and data 
processing than complete state tomography (see Meth- 
ods). 

We focus on a class of multipartite states routinely 
prepared in our experiment, the genuinely entangled n- 
qubit state \GHZ) n = i (|0)® n + |l)® n ) which yields 
the maximum value of our DIEW parameter 



T7nax 
nm 



2m"~ 1 cos(7r/2m) > B n 



when the following m measurements act on each qubit j 
in the x-y plane of the Bloch sphere 
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7T 

— for sj = 0, . . . , m — 1, (2) 



where a^} y ,z are the Pauli matrices for qubit j (see Sup- 
plementary Information for the optimization procedure) . 
For a given measurement s on all qubits, this measure- 
ment setting is then equivalent to rotating each qubit 

with exp(— i-^-ai^) for every j, followed by a collec- 
tive rotation of all qubits with exp(— i| J2j Vx^) where 
6 = 7r/2, and a measurement in the computational basis. 
Our witness is indeed device- independent: if these op- 
erations are not implemented as expected, because the 
(f) Sj take unprescribed values, or 8 ^ 7r/2, and conse- 
quently the measurement operators M 3 | s change, then 
by construction [3J the bound I nm < B nm still holds for 
all biseparable states whatever the actual measurement 
operators . are, thus avoiding false positive entan- 
glement detection. 

In the experiment each qubit was encoded in the in- 
ternal electronic Zeeman levels ASi/ 2 (m — —1/2) = |1) 
and 3D 5/2 (m = -1/2) = |0) of a 40 Ca+ ion. Three, 
four or six ions were confined to a string by a linear Paul 
trap and cooled to the ground state of the axial center- 
of-mass mode |25j . The entanglement among qubits was 
realized through a M0lmer-S0rensen |26l [27] entangling 
operation using a bichromatic light field collectively illu- 
minating all ions [T3] [e.g., Fig. la]. Single-qubit rota- 
tions exp(— i^-ai^) were driven by a far off-resonantly 
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detuned laser pulse focused on the target ion j and in- 
ducing an AC-Stark effect, where <j)j is determined by the 
detuning and the pulse duration [e.g., Fig. lb]. Collec- 
tive qubit rotations exp(— <Tx) were driven by a 
laser pulse exciting the qubit transition, where 9 is given 
by the Rabi frequency of the qubits and the duration 
of the pulse [e.g., Fig. lc]. At the end the qubits were 
measured in the computational basis through an electron- 
shelving technique scattering light on the dipole transi- 
tion iSi/2 O 4P]/ 2 . The scattered light was detected 
with a CCD camera that resolves each ion's fluorescence. 

In this setup, the subsystems may indeed leave the 
qubit subspace, albeit with an expected very small prob- 
ability, during initialization, preparation and measure- 
ment of the entangled state (see Methods). Using a 
DIEW allows us to discard these considerations safely. 

In practice we observe that applying a qubit rotation 
on an ion produces an unintended rotation on neigh- 
boring ions (see Fig. lb and Methods). Because this 
crosstalk effect is not included in the device-independent 
framework, which assumes all measurements to be de- 
fined locally, we analysed the impact that these crosstalks 
can have on the bound B nm . For simplicity, however, 
we only considered a crosstalk upper bound e. In our 
system we found that e = {1.6,5,4}% ± {0.1,1,1}% 
for {3, 4, 6} ions. We then computed numerically the 
maximum impact AI CT = Il isep — I^Zsep * nat crosstalks 
bounded by e can have on the biseparable bound for 
the settings intended for the experiment (CT stands for 
crosstalks) . Assumingthat the maximum contribution of 
the crosstalks to Eq. [yis given by AI CT , we update the 
bound to = B nm + AI CT . 

With three ions, we performed a DIEW measurement 
on a state prepared towards \GHZ) 3 as described above. 
For the m=2(3) settings DIEW of Eq. [2] each measure- 
ment set consists of 8 (18) measurements. To acquire 
significant statistics, we performed 6 (264) measurement 
sets where each measurement was done on 250 consecu- 
tive copies of the state (Table I) . For each set the order of 
the measurements was chosen randomly in order to avoid 
correlations between outcomes and unaccounted changes 
or drifts on the experimental conditions. 

To demonstrate genuine four- and six-partite entan- 
glement with DIEWs, we avoided a decay of coherence, 
quadratic in the number of qubits [TJ] (see Methods) , by 
preparing the decoherence-free GHZ states ^j(|0011) + 

1 1100)) and ^(|000111) + |111000>). These states have 
a coherence time of s» 300 ms, sufficient for our mea- 
surements. For a DIEW measurement of Eq. [2j the an- 
gles <j> 3j that maximize the violation for these states de- 
pend on the number of settings m and the number of 
qubits n (see Supplementary Information). For m = 2 
settings (sj = 0,1), a maximum violation is achieved 
when one half of the qubits is measured with an angle 
<f> Sj = SjT; and the other half with cf)' s , — + ~2 3 tt- 




Pulse duration (us) 



Figure 1. State and measurement characterization 
for the six-qubit DIEW measurement. a, Recon- 
structed density matrix p (absolute value) of six-qubit state 
~ 1000111} + e" 1 / 2 1 111000). The reconstructed state has a 
fidelity of 91.9(3)% with the target state, b, Single-qubit ex- 
citations while driving a single-qubit rotation of the fifth qubit 
(red markers) on a six-ion string, in this case, the crosstalk, 
due to imperfect focusing, is biggest for the fourth ion (black 
open circle markers) . The red line corresponds to a sinusoidal 
fit to the excitation of the fifth ion, and the black line to the 
fourth ion. The corresponding fits to the excitation of ions 
1, 2, 3 and 6 (markers dark blue, orange, green, and light 
blue) are not shown as the respective excitations are negligi- 
ble, c, Single-qubit excitations during a six-qubit collective 
rotation. All ions are homogeneously excited as the pulse du- 
ration is increased. The sinusoidal fit, shown as a black line, is 
done for all points. Error bars are ler where shown, otherwise 
smaller than the marker. 

For m = 3 settings (sj — 0, 1, 2), the corresponding an- 
gles are 4> Sj = + and cp' s , = ^^7r. The mea- 
sured DIEW parameters shown in Table I significantly 
violate the thresholds for genuine four- and six-partite 
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TABLE I. Summary of DIEW measurements with n = 3, 4, 6 qubits and m = 2, 3 settings^ 



n 


m 


B 










V = / cxp // max 


q (%) 


copics/mcas. 


meas./set 


sets 


copies 


3 


2 


4 


4.070(4) 


4.78(6) 


12 




0.84(1) 


15(1) 


250 


8 


6 


12k 




3 




10.542(9) 


12.39(1) 


136 


9V3 


0.7948(9) 


14.9(1) 


250 


18 


264 


1188k 


4 


2 


8 


8.43(8) 


10.42(6) 


20 


8^2 


0.92(1) 


19(1) 


50 


16 


60 


48k 




3 


18V5 


32.5(2) 


42.53(8) 


41 


27%/3 


0.909(7) 


23.7(6) 


50 


54 


60 


162k 


6 


2 


32 


34.0(5) 


40(1) 


5 


32^2 


0.88(3) 


15(3) 


5 


64 


21 


6720 




3 


162a/3 


294(3) 


374(3) 


19 


243%/3 


0.89(1) 


21(1) 


5 


486 


21 


51030 



a Errors in parenthesis, l<r, were derived from propagated statistics in the measured expectation values with 1000 Monte Carlo samples. 
B: biseparable quantum correlations bound, B CT : updated B to account for crosstalks, l ex P; measured DIEW, J max ; quantum DIEW 
bound, V: visibility, q = (I ex P - B CT )/I ex P: maximum admissible additional white noise. The magnitude of the violation, 
j-exp _ gCT ^ jg gj ven m un its of cr (gray column). 



TABLE II. Comparison of DIEW measurements with 
m = 2 and m = 3 settings for a similar number of 
copies of the state of n = 3, 4, and 6 qubitsQ 



n 


m 


jexp 


jexp _ b ct 


q (%) 


copies 


3 


2 


4.67(9) 


7 


13(2) 


4800 




3 


12.2(2) 


8 


14(2) 


4752 


4 


2 


10.5(2) 


10 


20(2) 


3840 




3 


41.7(6) 


14 


22(2) 


3240 


6 


2 


41(1) 


6 


17(3) 


1344 




3 


376(6) 


12 


22(2) 


1458 



a Errors in parenthesis, la, were derived from propagated 
statistics in the measured expectation values with 1000 Monte 
Carlo samples. B CT : updated biseparable quantum 
correlations bound to account for crosstalks, I ex P : measured 
DIEW, q = {I ex P - B CT )/I ex P: maximum admissible 
additional white noise. The magnitude of the violation, 
is given in units of a (gray column). 

entanglement. 

To compare the two- and three-settings DIEW mea- 
surements with the same number of copies, we chose 
random subsets of copies from each measurement of the 
above data. The results in Table II show a higher vio- 
lation for the three settings measurements with four and 
six ions. This indicates that the detection of entangle- 
ment through our witness is made easier by using more 
settings. Moreover, the maximum amount of white noise 
q that can be mixed with the observed statistics before 
loosing the violation indicates that more states are de- 
tected by the 3-settings than the 2-settings DIEW (see 
also Table I). 

Because DIEW inequalities with m = 2 settings are 
equivalent to n-partite Svetlichny inequalities [281 I29j . 
they detect not only multipartite entanglement but also 
multipartite quantum nonlocality [3J. However, for m > 
3 settings, the DIEW inequalities are not Svetlichny 
inequalities: Setlichny inequalities witness both, mul- 
tipartite entanglement and quantum nonlocality, but 
DIEW inequalities, in general, only witness multipar- 
tite entanglement (multipartite entanglement is neces- 
sary in order to show multipartite quantum nonlocal- 
ity) [3J. Multipartite entanglement can thus be detected 
device-independently even when no Svetlichny inequal- 



ity can be violated. This is indicated by the results 
of Table II for n = 6: For some amount of additional 
noise, the violation of the 2-settings inequality disap- 
pears while the 3-settings inequality can still be violated. 
Note that violation of a three-partite Svetlichny inequal- 
ity was previously demonstrated with photons in absence 
of crosstalks [3D], but leaving the locality and detection 
loopholes open. 

With today's technology, the different parts of an 
entangled system cannot always be well separated 
from each other. However, clear separation of the 
parts is certainly a desired feature in the long term, 
since it is necessary to take full advantage of quantum 
entanglement. As shown above, entanglement can be 
demonstrated in these circumstances without relying 
on a precise description of the measurement devices 
used, a knowledge which is never available completely. 
The achieved demonstration was thus more robust 
than conventional entanglement demonstrations and 
nevertheless simple. In order to guarantee the absence 
of crosstalks, we look forward to a demonstration with 
clearer separation between the parties or closing the 
locality loophole. On the theoretical side, further work is 
needed to conceive device-independent tests of resource 
states sufficient for interesting quantum simulation or 
communication protocols. 
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METHODS SUMMARY 

DIEW and tomography. For n parties and m possi- 
ble measurement settings per party, the total number of 
measurements settings required for the witness in Eq. 
is 2m n ~ 1 whereas tomography requires (d 2 — l) n measure- 
ments, where d is the dimension of each party or subsys- 
tem pp. For m = 3 and qubits (d = 2), for example, this 
witness requires two thirds the number of measurements 
required for full quantum state tomography (this advan- 
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tage becomes more evident as the dimension of the sys- 
tem increases) . In addition, the required post-processing 
is trivial in comparison with state reconstruction meth- 
ods. Naturally, considering m = 2 gives the minimum 
total number of measurements required, but as shown in 
the main text, considering m — 3 should result in more 
robust violations. For m — 4 and d — 2, however, the 
DIEW measurement requires more measurements than 
tomography. 

Experimental leak of qubit subspace. The qubits 
are initialized into the state 4S , 1 / 2 (rrij = —1/2) = |1) 
with a 99.9% probability, with the remaining population 
in the ASi/ 2 (mj = +1/2). Also, during the entangling 
operation, the motional ground state can be unintention- 
ally abandoned for higher motional states. Furthermore, 
our state detection protocol identifies as population in 
the state |1) any observation of fluorescence on the dipole 
transition, including ASi/2(mj — —1/2) and higher vi- 
brational states. 

Experimental crosstalks. In practice we observe 
that applying the rotation exp(— i^-ai^) on ion j pro- 
duces an unintended rotation ipk on qubit k because of 
residual light from imperfect focusing [25] (e.g., Fig. lb). 
This imperfection affects the net rotation applied on ion 
k, so that the measurement applied on ion k depends 
slightly on the one applied on ion j. These crosstalks 
can be quantified by the matrix Cjk = tfk/^Pj, and the 
operation intended only on ion j is then described by the 
operation on all ions exp(— i J2k ^Cjkcr^). For simplic- 
ity, we only considered a crosstalk upper bound e > Cjk- 

Collective decay of three-qubit state. Although 
the measured DIEW parameters shown in Table I for 
n = 3 parties do violate the threshold for genuine three- 
partite entanglement, the observed visibilities of about 
80% are incompatible with the 97% fidelity reported ear- 
lier Q3] . This disagreement is due to the collective decay 
of the state during the DIEW measurement. Although 
the single-qubit coherence time in this experiment was 
salO ms, the coherence time of the state \GHZ)% was only 
«2 ms. Because each single-qubit rotation took w 100 /is 
for a 2tt rotation (for collective rotations « 20 ^ls/2tt), 
and for each measurement the qubits were rotated se- 
quentially, the state was significantly decaying during the 
DIEW measurement. 
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I. EXPLICIT EXAMPLE FOR THE EVALUATION OF 
THE DIEW PARAMETER J nm 

Here we show an explicit example of the terms involved 
on the DIEW parameter 

Inm= {~l) S,m E S + 

s=0 (mod m) 

£ {-l^-W^Eg (1) 

s=l (mod m) 

where Eg — J2r(~~^Y P(r\s) is the n-partite correlator; 
s = Yl j Sj and r = . rj are the sums over all indices 
used in the measurement vector s and the outcome vector 
f (the colors will be evident in a moment). Obviously, 
for any settings s the following is true J2p P{r\s) = !• 

For n — 3 qubits and m = 3 measurement settings, 
there are 3 x 3 x 3 = 27 possible measurements as shown 
in the leftmost column in Table I. However, because the 
parameter I n= 3 iTn =3 involves only those settings for which 
s mod 3 is equal to or 1, only 18 of the 27 measurement 
settings are necessary. These entries, involving ,s modulo 
3 equal to or 1, are shown in red and green font colors 
in Table I. 

For our scenario where a measurement on each party 
has one of two possible outcomes, in our example with 
three parties, there are then eight possibilities, that is the 
vector of outcomes can take any of the following values: 

>re {000, 001, 010, 011, 100, 101, 110, 111}, 

where the vectors of possible outcomes are written here 
in boldface to distinguish them from the vector of mea- 
surement settings. 



For example, when the vector of measurement settings 
is s = (si, s 2 , s 3 ) = (0. 1.2), then ,s = 3 and s = mod 3, 
as shown in the row colored light gray in Table I. The 
3-partite correlator term E i2 = J2f(~ I) r -P(7i012) is 
therefore involved in the parameter ^33 as 

(_l)(o+i+2)/3£ ol2 = _[p(oo0|012) - P(001|012) 
+P(010|012) - P(011|012) 
+P(100|012) - P(101|012) 
+P(110|012) - P(lll|012)] 

In the scenario considered in which the measure- 



TABLE I Measurement settings considered in the evaluation 
of the DIEW parameter 133, see text. 
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merits are prescribed by angles, the settings are 

{^,^,^3} = {0 O ,01,0 2 }- 

Another example, when s = (si, 52,53) = (2,1,1), 
then s — 4 and s = 1 mod 3, as shown in the row col- 
ored dark gray in Table I. The 3-partite correlator term 
£-211 = S»-(~l) r ^(^l211) is therefore involved in the pa- 
rameter / 33 as 



(-1) 



( (2+i+i)-i)/3£ 2ii = _[p(ooO|211) - P(001|211) 
+P(010|211) -P(011|211) 
+P(100|211) -P(101|211) 
+P(110|211) - P(lll|211)] 

Again, in the scenario considered in which the mea- 
surements are prescribed by angles, the settings are 
{0 S1 ,0 S2 ,0 S3 } = {02, </>i, </>i }• For the 3-partite GHZ 
state, for example, the optimum angles calculated as de- 
scribed below, are 
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II. THE DIEW PARAMETER AS THE EXPECTATION 
VALUE OF THE QUANTUM OPERATOR I nm 

The value of I nm can be understood as the expectation 
value of the operator T nm on the measured quantum state 
p, that is I nm = Tr(I nm p). From the definition of the 
parameter I nm and the observed correlations P(f\s) = 



Tr 

as 



the operator can naturally be written 



2»m= E (-ir^Et-^i^ 

s=0 (mod m) r 



E (-!) (s - 1)/m < K 

s= 1 (mod m) r 



E E(- 1 ) r+s/m J 1 ^i Sj + 

(mod rn) r 

E Ei- 1 )^!^ 



s=l (mod m) r 

where M^., is the measurement operator acting locally 
on subsystem j corresponding to setting Sj and outcome 



numerical optimization. In general, this can be achieved 
by parametrizing the measurement settings M 3 r ^ s of 

each party j inside their respective Hilbert space W 
with a finite number of real variables. Gathering these 
variables for the measurement settings of all parties in 
a vector of parameters a then allows one to express 
the value of I nm as a function of a only, by letting 
P(f\s) = Tr(®] =1 M J r . ls Xa) p) in equation (1). The 

optimal settings are then given by M r , (a*) where 
a* = argmax I nm (a). 

In the case that the tripartite expression with m — 3 
settings J 3 3 is evaluated on a GHZ state 



\GHZ) = ^(|000) + |111)), 



for instance, the measurement operators corresponding 
to setting Sj of party j can be written as: 



ML = 1 + n,fJ , Ml, = 1 - ML , (2) 



where ||n|| = 1 and a — (o x , a y ,<T z ) is the vector of Pauli 
matrices. They can thus be parametrized by two real 
numbers . and ip{ . by letting: 



sin(tf^)sin(^.) 
cos(^) 



(3) 



18 parameters are thus enough to parametrize all mea- 
surements appearing in the ^33 expression. These can 
be gathered in the variable a — {$ r .\ g , V^..| s }s ,j — 
{^,^,^,i?i,...i?|,...,i?|,^,...,^l} 3 . 'The value of 
7 3 3 can then be expressed as a function of a di- 
rectly, by substituting P(r\s) = (GHZ\M^ lsi 0M^ S2 ® 
M^ s JGHZ) in equation (1). For instance, the prob- 
ability of observing outcomes 000 when all parties use 
their first setting is given by 



P(000|000) = (GHZ\Ml <g> Ml ® Ml \GHZ) (4) 



III. OPTIMAL MEASUREMENT SETTINGS FOR THE 
DIEW PARAMETER I nm 

For every quantum state p, the measurement settings 
that yield the largest value of I nm can be determined by 



After optimization of J33 over a, one finds that all 
values of i) 3 s . take the value tt/2. This means that the 
settings yielding the largest value of 13,3 upon measure- 
ment of a GHZ state lie in the x-y plane of the Bloch 
sphere. 



